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1 Introduction 



In a recent publication [1] we have suggested a relativistic field theory model of the 
deuteron. This model, being some kind of a cr-model [2], is based on the assumption 
that the physical deuteron state should be produced due to integration over low-energy 
proton-neutron fluctuations at energies restricted by the scale Ad. For simplification 
we have suggested a one-nucleon loop approximation for the integration over proton- 
neutron fiuctuations. Unfortunately, this approximation does not have a perturbative 
parameter and should be considered as effective one, like that suggested by Nambu and 
Jona-Lasinio [3]. In this case the scale Ad has the meaning of the cut-off. We have defined 
Ad in terms of the effective radius of the deuteron Td, i.e. Ad = 1/td- For the estimate 
of Ad we have applied the non-relativistic formula: td = (^d -^n)" ''^''^ — 4.319 fm 
[4,5], where £d — 2.225 MeV is the binding energy of the physical deuteron [5] and 
Mn = 938 MeV is the mass of the nuclcon. We used equal masses for the proton and 
neutron, i.e. Mp = Mn = Mn = 938 MeV. This corresponds to the chiral limit when 
masses of current u- and d-quarks vanish, i.e. mou = ^od = 0. Our estimate of Ad 
gives: Ad = l/ro = 46 MeV. 

The interactions of the deuteron with proton and neutron were described in terms 
of two coupling constants gy and qt, defining the interactions via vector and tensor 
nucleon currents, respectively. In the one-nucleon loop approximation we calculated the 
anomalous magnetic dipole moment of the deuteron kd in units of the nucleon magneton 
A*N = e/2M^ , where e is the proton charge, and Qd the electric quadrupole moment. 
Imposing the constraint kd — 0, being valid in the lowest approximation, we found the 
correlation between coupling constants gy and qt, i.e. Qt/qy — — •\/3/8. Then the 
coupling constant gy has been fixed in terms of the electric quadrupole moment Qd- 
Eventually we calculated in one-nucleon loop approximation the binding energy of the 
physical deuteron in terms of gy and Ad. The theoretical result obtained was in good 
agreement with the experimental value: (£D)exp — 2.225 MeV [5]. 

The physical nature of the coupling constants gy and g^ is connected with one-meson 
exchange. Of course, the vr-meson exchange should give the main contribution. This 
assumption has been confirmed by the magnitude of gy, i.e. gy ~ g'TrNN , where g'TrNN is 
the coupling constant of the TrNN-interaction. 

The application of the model to the calculation of processes of low-energy interactions 
of the deuteron encounters the problem of the unambiguous computation of one-nucleon 
loop digrams that should describe the interactions of the deuteron with other particles in 
the suggested model. 

Indeed, since these are fermion loops, most of them, contributing to processes like 
radiative neutron-proton capture n + p ^ D + 7, low-energy proton-proton reaction 
p + p— >D + e'^ + i/e and so on, depend explicitely on the shift of virtual momenta 
of nuclcons in the loop. This introduces substantial ambiguities interfering with the 
direct application of the model. However, free parameters appearing due to shifts of 
virtual momenta in nucleon loops must not be considered as free parameters of the model 
introduced especially for the computation of every nucleon diagram. 

In this paper we revise our model [1] and apply it to the calculation of the cross 
sections of the radiative neutron-proton capture n + p — D + 7 and two-proton fusion 
p + p— ^D + e^ + i^e. In order to compute unambigously the effective anomalous 
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magnetic and electric quadrupole moments of the dcutcron we impose constraints on 
the computation of one-nucleon loops, caused by the requirement of electromagnetic 
gauge invariance of contributions for individual loops. This distinguishes our effective 
model from a field theory like QED, where the electromagnetic gauge invariance should 
be required for the complete set of diagrams in fixed order of perturbation theory. The 
requirement of electromagnetic gauge invariance applied to the individual nucleon loops 
allows one to fix ambiguities and use one-nucleon loops as well-defined quantum-field- 
theory objects. In order to remove ambigutities of the computation of the one-nucleon 
loop diagrams describing the effective Lagrangian of the radiative neutron-proton capture 
we turn to the selection rules. In the case of the computation of the effective Lagrangian 
responsible for the low-energy two-proton fusion p + p— s>D + e+ + z/eWe can remove 
ambigutities of the computation of the one-nucleon loop diagrams by using the gauge 
invariance under gauge transformations of the deuteron field. 

The present paper is organized as follows. In Sect. 1 we adduce the starting assump- 
tions that we have put in the foundation of the model, and phcnomcnological parameters 
in terms of the parameters characterizing the physical deuteron. In Sects. 2 and 3 we 
give a detailed derivation of the effective Corben-Schwinger and Aronson Lagrangians, 
respectively, describing electromagnetic interactions of the deuteron in terms the magnetic 
and electric quadrupole moments. In Sect. 4 and Sect. 5 we apply the relativistic field 
theory model of the deuteron to the computation of the amplitudes and cross sections of 
the radiative neutron-proton capture n + p — D + 7 and the low-energy proton-proton 
fusion p + p— ^D + e^ + i^e, respectively. In the Conclusion we discuss the obtained 
results. 

2 The deuteron structure in the one-nucleon loop ap- 
proximation 

In Ref. [1] we have suggested a relativistic field theory model for the deuteron as a 
bound state of proton and neutron. The main idea, that has been put into the foundation 
of the approach, has been adopted from the Nambu-Jona-Lasinio (NJL) model [3]. As 
has been suggested by Nambu and Jona-Lasinio [3], the physical deuteron should appear 
as the proton-neutron collective excitation due to a local four-nucleon interaction and 
vacuum fluctuations of nucleons, calculated in one-nucleon loop approaximation, produc- 
ing both a kinetic term of the physical deuteron and its interactions with other particles. 
Unfortunately, the direct realization of the Nambu-Jona-Lasinio program, i.e. the use 
of the local four-nucleon interaction, would lead to a strongly bound localized proton- 
neutron state, whereas a physical deuteron is weakly bound and rather "extended" object. 
Therefore, we have borrowed from the NJL model only the admission of the applicability 
of the one-loop approximation. This is used for the computation of observed parameters 
for physical states in leading order in the long-wavelength expansion [6-9]. 

In our description of the physical deuteron, being some kind of the cr-model, we 
use the following scheme: We start from the Lagrangian of an unphysical deuteron field 
D^^^x), considered as a bound proton-neutron state at zero binding energy and with a 
mass equal to the sum of the proton and neutron masses. Then in the one-nucleon loop 
approximation and leading order in the long-wavelength expansion we obtain an effective 
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Lagrangian of a physical deuteron field describing a physical deuteron with observable 
binding energy (£D)exp = 2.225 MeV [5], anomalous magnetic dipole moment (fi;D)exp = 
—0.023, determined as = A^d — /^p — /^n, where /id = 0.857, /Xp = 1 + Kp = 2.793 and 
/^n = ft^n = —1.913 are the magnetic dipole moments of the physical deuteron, proton 
and neutron, respectively, and the electric quadrupole moment (QD)exp = 0.286 fm^ [5]. 
The magnetic dipole moment of the deuteron is measured in nuclear magnetons /in = 
e/(2MN), where e and Mn are the electric charge of the proton and the mass of the proton 
and neutron, then and k,^ are the anomalous magnetic dipole moments of the proton 
and neutron. In the case of the neutron the total magnetic dipole moment coincides 
with the anomalous one. Below, to simplify matter, we neglect the proton-neutron mass 
difference and use equal masses of the proton and neutron, i.e. Mp = Mn = Mn = 
938 MeV. This should correspond to the chiral limit approximation with zero masses of 
current u- and d-quarks, i.e. mon = mod = 0. 

In this connection we would like to refer to the paper written by Sakita and Goebel [10], 
where, for the aim of the investigation of the low-energy limit of the photodisintegration, 
the deuteron has been described in terms of a local field operator. 

The Lagrangian of the unphysical deuteron field D^^\x), which interacts strongly with 
the proton p{x) and neutron n{x) fields, reads 

-igy [p{x)'y'^n''{x) - n{x)'y'^p''{x)] Dj^\x) - 

-igy[p''{x)j''n{x) - n''{x)-f^'p{x)]Dl^°\x) + (2.1) 

Mq 

+ ^[pc(x)(j^"'n{x) - n'ix)(j^"'p{x)]Djj^\x) + 
Mo 

+p(x)(i7^9^ - Mn)p(x) + n{x){i-f''d^ - Mt^)n{x). 

Here D^^^{x) = d^D'^\x) — d^D^^\x), Mq = 2Mn is the mass of the unphysical deuteron, 
then ip'^{x) = Cip'^{x) and ip'^^x) = ip'^{x)C, and a'^" = 7*^] ■ The operator 

C denotes charge conjugation, T is a transposition, gy and g^ are the phenomenological 
constants which will be fixed below. We assume that the coupling constants gy and 
gT are caused by one-meson exchanges. The vr-meson exchange should give the main 
contribution. 

The guiding principle for the construction of the Lagrangian ( p.l|) has been only the 
fact that the deuteron is a bound state of the proton and neutron with spin one [2]. 

In order to compute the electric quadrupole and anomalous magnetic dipole moments 
of the physical deuteron we have to include the interaction the unphysical deuteron field 
Djj^\x) with the electromagnetic field. Having performed this inclusion by a minimal 
coupling we obtain the Lagrangian 

Ctot{x) = C,t{x) + Cei{x), (2.2) 

where 

C4x) = -ieDlf{x)A^{x)D^''^''{x) + ieDfXx)A^{x)D^^''^''{x) + 
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c 



{x)cT^''p'{x)]A^{x)Dl^\x) 



{x)a^''p{x)]A^{x)Dl^^\x) 



(2.3) 



A 



Dt(o)(^)Z^(o) 



{x)Fa^{x) - ep{x)'y^p{x)Af,{x) + 





p{x) (t'^''p{x) F^^{x) + ie 



■n 



n(x) a>"'n{x) F^^{x) + 0{e^). 



Here Fa^{x) = daA^{x) — d'^Aai^x) and Afj,{x) are the electromagnetic field strength 
tensor and the electromagnetic potential, respectively. Also we have added the Aronson 
interaction [11] describing the electric quadrupole and anomalous magnetic moments of 
the unphysical deuteron: = 2A/Mq and k^co) = A [11]. 

Now let us give some arguments for the justification of the application of the Aron- 
son interaction. Recall, that the electric quadrupole and the anomalous magnetic dipole 
moments characterise the inhomogeneity of the structure of the unphysical deuteron rep- 
resented in terms of the field Dj^^^x). The availability of such an inhomogeneity is not 
obvious beforehand and does not follow from the Lagrangian (|2.1| ) via a minimal coupling 
inclusion of the interaction of the unphysical deuteron with an external electromagnetic 
field. Therefore, the inclusion of nonzero Qdw and Kj^io) is an additional assumption 
concerning the internal structure of the unphysical deuteron field. Since, the minimal 
coupling leads to the appearance of electromagnetic interactions linear in Dj^J(x), the 
required interaction of the unphysical deuteron field with an external electromagnetic 
field, caused by electric quadrupole and anomalous magnetic dipole moments, should be 
of order 0{D'jj'J{x)D^^p{x)). This does not contradict to the main principles of Classical 
Electrodynamics [12] and gives the required interaction in the form, being irreducible 
with respect to those interactions induced by a minimal coupling. The Aronson effective 
Lagrangian provides a simplest form of an interaction of order 0{Dj^J{x)Dl^^{x)) [11]. 

The inclusion of the Aronson interaction can be also justified by the assumption that 
the electromagnetic interactions induced by electric quadrupole and anomalous magnetic 
dipole moments of the unphysical deuteron field with an external electromagnetic field 
become available only in next-to-leading order in long-wavelength expansion, i.e. in 
powers of gradients of the deuteron field Djj^J(x). In leading order in the long-wavelength 
expansion, linear in Dj^J{x), the interaction of the unphysical deuteron with an exter- 
nal electromagnetic field, caused by electric quadrupole and anomalous magnetic dipole 
moments, vanishes. 

We have taken into account the anomalous magnetic dipole moments of the proton 
and neutron that play an important role for the radiative neutron-proton capture n + p 
— >^ D + 7. The inclusion of A and the anomalous magnetic dipole moments of the proton 
and neutron distinguishes the present model from that given in [1] . 

Some of the results and assertions adduced in this section repeat those given in [1]. 
This concerns mainly the computation of the binding energy sd and the discussion of 
the validity of the approximation having been applied to the computation of one-nucleon 
loop diagrams. This repetition is caused by the consideration of the completeness of the 
exposition, i.e. in order not to relegate a reader to [1] for consultation, and the change 
of the magnitudes of phenomenological parameters of the model. The former is due to 
different procedures of the removal of ambiguities of one-nucleon loop diagrams, induced 
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by shifts of virtual nucleon momenta, we use here and we have used in Ref. [1]. In the 
present version all shift ambiguities are fully fixed by requirements of gauge invariance. 

Following [1] we compute the effective Lagrangian of the physical deuteron field in 
the one-nucleon loop approximation (see also [6-9]). The nucleon diagrams describing 
the contributions to the kinetic term of the deuteron field and leading to the appearance 
of the non-zero value of the binding energy are depicted in Fig. 1. By calculating these 
diagrams in leading order in the long-wavelength expansion we get 



The ultra-violet cut-off A restricts the 3-momenta of fluctuations of virtual nucleons 
taking part in the formation of the physical deuteron field. One should expect that A is 
connected with the region of localization of a wave packet procreated by fluctuations of 
virtual nucleons, i.e. Ar • A 1. We shall specify the value of A below. 

Now we discuss the problem of the applicability of the long-wavelength expansion for 
the calculation of the one-nucleon loop diagrams in Fig. 1. It is well-known that a two- 
body S'-wave bound state being denoted as D with a reduced mass and a binding 
energy Ed is localized in the region restricted by rn = I/a/^d^n [4]. This quantity can 
be considered as the effective radius of the bound state [5]. In the case of the physical 
deuteron we have = 4.319 fm [5]. This value exceeds three times the effective radius 
of nuclear forces tn = = 1.462 fm, where M^^ = 134.976 MeV is the pion mass [13]. 
Hence the physical deuteron looks like a rather "extended" bound state. It should be 
obvious that such an "extended" bound state can be formed at the expense of the main 
contributions of long-wavelength fluctuations of the bound proton and neutron. It implies 
that we can expect a cut-off A satisfying the inequality A ^ Mn. It should be obvious 
that A is to be identified with Ad = l/ro = 45.688 MeV, i.e. A = Ad = 45.688 MeV. 

The infinitesimality of the derivatives d^D^^"* (x) necessary for the validity of the long- 
wavelength expansion can be justified as follows. We assume that the deuteron field 
D^^\x) is a bound state of a proton and a neutron at zero binding energy. As has 
been mentioned above, this means that the field Df^'^ (x) is localized in the region whose 
upper boundary rD(o) goes to infinity. Obviously it results in a smooth variation of the 
field D'^lf'{x) at the scale of the effective radius of the physical deuteron. Thus we have 
adduced some arguments on behalf of the validity of the long-wavelength expansion which 
has been applied for the calculation of the effective Lagrangian ( p.4|) . 

Now we proceed to the computation of the binding energy of the physical deuteron. 




(2.4) 



where Ji(Mn) 



and J2(Mn) are the following divergent integrals 




(2.5) 
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First, following [1], we introduce the field of the physical deuteron 



D,{x) = Z'^/'D^^\x), (2.6) 

where 

is the wave-function normalization constant. After the renormalization (7) we get the 
effective Lagrangian of the free physical deuteron field D^{x) [1] 

4S(x) = -]^Dl^{x)D^'^{x) + MlDl{x)D^{x), (2.8) 

where Mq = Mq — sd is the mass of the physical deuteron and is the binding en- 
ergy. In the one-nucleon loop approximation and in leading order of the long-wavelength 
expansion the binding energy Ed is determined by the expression 

= ^ Ji(M,) + - 4^v^T + 2,1 j^^^^^^ (2.9) 

D TT^ ik/M 27r^ 



For the derivation of Eq. ( p.9|) we have used the inequality 

9w - Q9y9T + 3,T 



3 7r2 



J2(Mn) < 1 (2.10) 



being consistent with the inequality A^ ^ Mn discussed above. For Ad <^ Mn Eq- 
reads 



r - A3 1 
— 



9 7r2 M2 



12 ^ 6 /^y 

5 Uv/ 5 [gyl 



(2.11) 



It is seen that the binding energy of the physical deuteron is expressed in terms of three 
phenomenological parameters of our model: Ad, Qy and qt- Then we have to fix the 
coupling constants gy and to obtain the theoretical value of the binding energy. 

Eqs. and (|2.11|) have been calculated in [1] . We need to reproduce them here for 



two reasons. First, we give below the new relation between the coupling constant gy and 
the electric quadrupole moment Q^. Second, in the Conclusion we shall discuss Eq. ( p.ll| ) 
setting (7t = 0. For these reasons it is convenient to have Eqs. ( ^.9| ) and ( |2.11| ) at hand 
and not to relegate a reader to Ref . [1] . 

In order to fix the coupling constants gy and (7t "we suggest to calculate the anomalous 
magnetic dipole kd and electric quadrupole Qd moments in terms of gy and We admit 
that Kj) and Qb appear in the one-nucleon loop approximation, and only at the expense 
of interactions taken into account in the Lagrangian (|2.2| ) [1]. The complete set of one- 
nucleon loop diagrams is depicted in Fig. 2. The non-trivial contributions come from the 
diagrams in Figs. 2c, d. The detailed computation of the diagrams depicted in Figs. 2c, d 
is given in Sects. 2 and 3, respectively. Here we only adduce the results. 
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The effective Lagrangian, defined by the one-nucleon loop diagrams in Fig. 2 and de- 
scribing the electromagnetic interactions of the physical deuteron field via the anomalous 
magnetic and electric quadrupole moments, is given by 



^e^Dlix)D,ix)F^'^ix) + 



+ ie< 



2g^ 
37r2 



9 

1 - -f/tr 



1 - 



-A 



(2.12) 



The Lagrangian ( p. 12 ) contains only physical deuteron fields. The first term in £^g(x) 
is the Corben-Schwinger interaction [14], while the second term presents the interaction 
which has first been introduced by Aronson [11]. These interactions describe the anoma- 
lous magnetic dipole and electric quadrupole moments of the charged vector field. It 
should be emphasized that we have neglected the divergent contributions which are small 
in comparison with the convergent ones due to the restriction ( 2.10| ). 

Note that the effective Lagrangian (|2.12|) differs from that obtained in Ref. [1] by 
the additional contributions caused by anomalous magnatic moments of the proton and 
neutron, and the phenomenological coupling A induced by the Aronson interaction. If 
one would set = A = 0, one would find a discrepancy in the main terms. The 

former is due to different procedures of the removal of ambiguities of one-nucleon loop 
diagrams, induced by shifts of virtual nucleon momenta, we use here and we have used 
in Ref. [1]. In the present version all shift ambiguities are fully fixed by requirements of 
gauge invariance. 

The anomalous magnetic dipole moment kd measured in units of the nuclear magneton 
/^N = e/(2MN) and the electric quadrupole moment Qd are given by [10,13] 



kd = - 



4 7r2 



37r2 



9 

1 ■^■('^P '^n 




A 



(2.13) 



37r2 



9 

1 - -f^r 



A 



M2 



(2.14) 



The experimental value (ftD)exp = —0.023 of the anomalous magnetic dipole moment 
of the deuteron is small compared with the magnetic dipole moment of the deuteron 
fj,D = 0.857. Therefore, we suggest to set = giving 



27r2 



A - 



37r2 



1 - 



9 



1 - 



9^T, 



(2.15) 



In this case the observed magnitude (KD)exp = —0.023 can be the matter of next-to- 
leading order corrections to the accepted approximation, for example, chiral corrections, 
etc. 

By applying Eq. ( |2.15| ) we can express the electric quadrupole moment Qd in terms 
of the coupling constant gy only, i.e. 



Qd 



Ma 

7r2 



(2.16) 



8 



The relation between and Qd, represented by Eq. ( p.l6|) , differs by a factor 2/3 from 
tfiat obtained in Ref. [1]. 

Such a discrepancy has the following explanation. The relation Eq. ( p.l6| ) is caused by 
effective couplings of the deuteron to electromagnetic field, i.e. the Corben-Schwinger and 
Aronson interactions. These interactions are defined by one-nucleon loop diagrams the 
computation of which depends substantially on the shift of virual nucleon momenta (see 
Sects. 3 and 4). The mathematical procedure of the removal of such ambiguities which 
has been used in Ref.[l], unfortunately, has left room for residual ambiguities. In the 
present version ambiguities induced by shifts of virtual nucleon momenta are fully fixed 
by requirements of gauge invariance under gauge transformations of the electromagnetic 
and the deuteron fields. The former has led to the factor 2/3 in the relation between Qd 
and g'^. 

By using the experimental value (Qd)cxp = 0.286 fm^ and Md — 1876 MeV we can 
estimate the value of g^- This gives: gy = ±11.3. Without loss of generality we can use 
the positive sign, i.e. 

gy = n.3. (2.17) 

The coupling constant gy satisfies the relation gy ~ (?,rNN , where gn^'N = 13.4 ±0.1 [4] 
is the coupling constant of the vrNN-interactions. Thus the magnitude of gy corroborates 
our assumption that the phenomenological interactions of the deuteron and the proton 
and neuteron, given in the Lagrangian ( |2.1| ), are caused by one-meson exchange and 
the TT-meson exchange gives the main contribution. This admission can be justified by 
comparing the effective radii of one-meson exchanges. The radii of pion, a (660) -meson 
[6,9] and p (770) -meson exchanges are defined in terms of their masses, i.e. r^^ = l/M^^ = 
1.462 fm at = 134.976 MeV [13], = = 0.30 fm at = 660 MeV [5,9] 

and rp = l/Mp = 0.27 fm at Mp = 770 MeV, respectively. The radii of the a (660) 
and p (770) -meson exchanges are much smaller than the radius of the vr-meson exchange. 
Therefore, these interactions of the proton and neutron with a (660) and p (770) mesons 
should contribute perturbatively to the deuteron problem. 

We have put = as the experimental value (KD)exp = — 0.023 is small enough in 
comparison with the magnetic dipole moment of the deuteron /id = 0.857. We assume 
that the nonzero value of kd can be obtained perturbatively by taking into account nonzero 
values of current quark masses, for example, within Chiral perturbation theory at the 
quark level (CHPT)g [9], based on the extended Nambu-Jona-Lasinio model with linear 
realization of chiral U{3) x U{3) symmetry, and interactions of the proton and neutron 
with a (660) and p (770) -mesons within the one-meson exchange approximation. 

From Eqs. ( |2.11| ) and ( |2.16| ) we can express the binding energy eu in terms of the 
electric quadrupole moment Qd and the ratio gy/gT 



— — VD Aj-) 



^ _ 12 fgA _^ 6 /gr 
5 \gvl 5 Uv, 



2- 



(2.18) 



This formula agrees with the statement that the main contribution to the binding energy of 
the deuteron comes from the tensor forces producing the non-zero value of the quadrupole 
moment Qd [3]. For the computation of the magnitude of the binding energy we have to 
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fix the ratio Qt/ g\ that is still free. Putting Qt/ g\ 
of the binding energy 

= 2.273 MeV. 



— w3/8 [1] we obtain the best fit 



(2.19) 



The accuracy of the fit makes up 2%. Thus, we have fixed all phenomenological parameters 
and described all quantities, characterizing the physical deuteron. 

The total effective Lagrangian of the physical deuteron describing strong and electro- 
magnetic interactions of the deuteron reads 



]^Dl,{x)D^^{x) + MlDl{x)D^{x) - teDl{x)A^{x)D'' {x) + 



+ieD^^{x)A^'{x)D^''{x) +ie 
1 



27r2 



DUx)D,{x)F^'''{x) 



-le 



9\ 



27r2 



Ml 



D (x)D'''' x)R/ x) 



-ig\\p{x)'^^rf{x) — n{x)'y'^p'^{x)]D^{x) — 
-igy[p^{x)'y^n{x) — n'^{x)^'^p{x)]DUx) + 



(2.20) 



9t r- 



Mr 



\p{x)a^''n''{x) - n{x)a^''p''{x)]D^^{x) + 
[p''{x)(T^''n{x) - n''{x)a^''p{x)]Dl^{x) + 



+ie'^[pix)a^"'n''{x) - n{x)a>'''p%x)]A^{x)D,{x) - 

-ie'^[p^{x)a>"'nix) - n-(x)a'^Xx)]A^(x)Z^i(x) + 0{e^) + .... 
Mr, 

The ellipses stand for interactions of the proton and neutron with other fields such as 



photon, pions, etc. The effective Lagrangian (|2.20|) differs from that obtained in Ref. [1] 
by the coupling constants of the effective interactions of the physical deuteron with the 
electromagnetic field describing the anomalous magnetic dipole and electric quadrupole 
moment Qj). 

We have to underline that the suggested field theory model of the deuteron is applicable 
only at the low-energy limit. Thereby all interactions of the deuteron with other hadrons 
should run via one-nucleon loop exchange. This is due to the one-nucleon loop origin of 
the deuteron in our model. This assertion is very similar to the approximation accepted 
within the NJL model, where all interactions of hadrons run via one-constituent quark- 
loop exchange [6-9]. However, if the one-constituent quark loop approximation in the 
extended NJL quark model [9] can be justfied in large N (number of quark colours) 
expansion, where the perturbation theory is developed in powers of 1/N. In the case 
of the relativistic field theory model of the deuteron for the justification of the one- 
nucleon loop approximation we do not have any perturbative parameter. Therefore, our 
approximation cannot be justified like that in the extended NJL quark model and should 
be considered as an effective one. For the confirmation of the validity of this effective 
approximation we can only refer to Ref. [1], where the two-nucleon loop contributions to 
the binding energy have been calculated and found much less than the one-nucleon loop 
contributions. 
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We deem that most likely the deuteron cannot be inserted in an intermediate state of 
any process of low-energy interactions. This is connected with a very sensitive structure 
of the deuteron as an "extended" bound state with a small binding energy. The repre- 
sentation of such a state in terms of any local quantum field is rather limited. The latter 
entails an undetermined character of the description of the deuteron in intermediate states 
in terms of Green functions of these local fields. 

Now we can proceed to the computation of the effective interactions describing the 
anomalous magnetic and electric quadrupole moments of the deuteron. The complete 
set of one-nucleon loop diagrams describing effective interaction of the deuteron with 
the electromagnetic field is depicted in Fig. 2. In leading order in the long-wavelength 
expansion the diagrams in Figs. 2a,b do not contribute to the anomalous magnetic and 
electric quadrupole moments. The contributions of these diagrams are divergent and can 
be removed by the renormalization of the electric charge of the deuteron. The nontrivial 
contributions to the effective anomalous magnetic and electric quadrupole moments are 
defined only by the diagrams in Figs. 2c, d. 



3 The effective Corben— Schwinger interaction 

In this Section we give the detailed derivation of the effective Corben-Schwinger 
Lagrangian defined by the one-nucleon loop diagram in Fig. 2c and describing the effective 
electromagnetic interactions of the deuteron. The effective Lagrangian of the diagram in 
Fig. 2c is defined 

J d^x£Fig.2c(x) = J d^x J ^ ^ {27!)^^ ^^^^^ -Pl(^i) -4^(3^2) X 

2 

X e-^^^-^i e-^'^^-^^ ^^(k, + k2)■x ^ k2\ Q) , (3.1) 

4 TT^ 

where 

J^''^{Kk2;Q) = (3.2) 

d^k { g I ^ 1 „ 1 1 
— tr<^ 7^ . ^ 7" . . ^ 7^ . . . ^ . 

^^^^ I M^-k-Q - k - Q - ki - k - Q - ki - k2 } 

The 4-vector Q = aki + bk2, where a and b are arbitrary constants, displays the 



dependence of the k integral in (|3.1| ) on the shift of the virtual momentum. This ambiguity 
of the computation of the integral over k, which has been found by Gertsein and Jackiw 
[15], is used to remove undesirable contributions and make the effective Lagrangian gauge 
invariant. We use the fields of the physical deuteron. This is because the renormalization 
( p.6|) introduces divergent terms that are small compared with the convergent ones we are 
following for the computation of the Corben-Schwinger Lagrangian. 

In order to display the Gertsein- Jackiw ambiguity we follow the Gertsein- Jackiw 
method and compute the difference 

Sj''''^{k,,k2;Q) = J^'^^ik^^k^-.Q) - J^"^(fci, ^2; 0) . (3.3) 
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In accordance with the Gertsein-Jackiw method the difference (3^) can be represented 
by the integral 



dx 



TT^i d k^ 



1 



tr<^ Y ^ ^ 7" 

[ M^-k-xQ M^-k-xQ -k] 



X 



Mn — k — xQ — ki — k2 . 



(3.4) 



This shows that the Gertsein-Jaclciw ambiguity is just the surface term. Following Gert- 
sein and Jackiw [15] and evaluating the integral over k symmetrically, we obtain 



5J^''^'{ki,k2;Q) = -2 [' dx lim (^-^tT{^P (M^ + R + xQ)^ x 

Jo R^oo \ 



x{M^ + R + xQ + h) 7" {M^ + R + xQ + ki + h)}) ■ 



(3.5) 



The brackets < . . . > denote the averaging over directions of the 4-vector R. Due to the 
limit R ^ 00 we can neglect all momenta with respect to R. 



SJ^''''ik,,k2;Q) 



2 lim l^tTi^^RrRYR} 



(3.6) 



Averaging over directions of the 4-vector R 



we obtain 



lim 

ij— »oo 



R^R^R'^RP 



1 

24 



1 



(^A^^.P + g'^g''" + g^'g 



Xp „'P^\ 



{Q^'gfP + Q'^g^'^ + Q^'g 



,13 fia 



(3.7) 



(3.^ 



Thus the surface ambiguity noticed by Gertsein and Jackiw contains only finite contribu- 
tions. 

Now we should proceed to the computation of J^"^{ki, Q)- In order to pick up the 
ambiguity connected with Q one cannot apply the Feynman method of the computation 
of momentum integrals like ( p.2| ). This method involves the merger of the factors in the 
denominator with the subsequent shift of virtual momentum. On this way one can lose 
the Q-depenedence that is due to the shift at the intermediate stage. Thereby we have 
to compute the integral over k without any intermediate shift. One can carry out this by 
applying a long-wavelength expansion and keeping to the leading terms. 



d^k 



7tr<^7 



Mn + + g 

M2 - fc2 



1 + 



2k-Q 



7^ 



M^ + k + Q + ki 
- k^ 



X 
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X 



1 + 



2k-{Q + ki] 



^Mt^ + k + Q + ki + k2 



^ Ml-k^ 



1 + 



2k-{Q + ki + k2) 



Ml-k^ 



d^k 



TlH (M2 - A; 2) 3 



tr {M^7^(^ + g)7''7'' + M27^7"(^ + Q + ^1)7'' + 



+m27V7''(^- + Q + k + k) + 7''(^ + Qh'^ik + Q + ki)Y{k + Q + k + k)} 
2k- {3Q + 2ki + k2 " 



1 + 



M^-k^ 



(3.9) 



- 1 
~ 2J 



+ 



(M2-fc2)2 (M2-A;2)3 

+7V7^(Q + fci + A;2)} + 

d^A; A; • (3g + 2A;i + k2) 



tr {7^Q7"7'' + 7^7" (<9 + ^1)7'' + 



tr {M^(7^A;7"7^ + ^'^^''kY + ^f^Tl^k) + 



' (M2-A;2)4 
V^7"^7^^} = J(?"''(A;i,A;2;g) + J^^^^{k^M;Q)- 

For the computation of J^^^{ki^k2]Q) it is sufficient to calculate the trace of the Dirac 
matrices and integrate over k 



J^^^^^{k^M\Q) = [1 + 2 J2(Mn)] [(Q + 2A;i + A;2)V' + 
+ (g + 2A;i + A;2)V" + (Q + 2A;i + k2Yg"^ - 
-2{ki + k2Tg^^ - 2k^,gn, 



(3.10) 



where J2{M^) describes a divergent contribution depending on the cut-off Ad. Due to the 
inequality Ad ^ M-^ we can neglect J2(Mn) with respect to the convergent contribution. 

7-/3"/ 
'(2) 

gives 



In order to compute J{-^^{k\^ k2] Q) we have to integrate ffist over k directions. This 



Jll'^^{k^,k2■,Q) = 



rd'^k 
J n^i 



{3Q + ki + k2)xx 



2(M2-A;2)4 6(M2-A;2)4 

Xtr (7^7'^7"7'' + ryl^ryOC^X^H ^/S^U^H^X^ _ 

hi + 6 J2(Mn)] [(3Q + 2A;i + A;2)V + (3g + 2A;i + A;2)V" + 

y 



(3.11) 



+ (3g + 2A;i + A;2)^^"'']. 
Here we have used the integrals 

d^k 1 



r (Tk. 
J n^i 



2M2 ' 



6M4 



(M2 - A;2)3 

f '^^^ ^ 
J ^ (M2 - A;2)4 

Summarizing the contributions, we get 

J^'^'ik^. k2; Q)^l {Q^g^" + Q^g''' + Q'^g'^^) + ^ [1 + ^ UMn)\ x 



(3.12) 
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x[{2h + A;2)V^ + {2h + hfg^''' + {2k^ + hYg'^"] + (3.13) 
+ [1 + 2J2(Mn)][-2(A;i + k^Yg^^^ - 2A;Vl- 

It is seen that the Q-dependence coincides with that obtained by the Gertsein-Jackiw 
method (|3.8| ). Due to the arbitrariness of Q we can arbsorb by the Q-term the terms 
having the same Lorentz structure. This brings up the r.h.s. of ( |3.13| ) to the form 

+ [- 2 {ki + k^Tg"^ - 2 fcf (7^"] . (3.14) 

Also we have neglected the divergent contribution. This approximation is valid due to 
the inequality Ad -C Mn. 

The effective Lagrangian £Fig.2c(a^) defined by (p.l4|) reads 

/:Fig.2c(x) = ^e||j[(3 - a)d''Di{x)Dp{x)A^{x) - 

-(3 - a)Dl{x)d^Dp{x)A''{x) - h DI{x) D p{x) d'' {x) - 

-{h - a) Di{x) Dp{x) A"{x) - (3.15) 

-(a - b)d'^ Di{x)D''{x)A(^{x) + 6D^(a;) 9^ D"(x) A^(x) + 

+ 3Dl{x)Dp{x) (d^A^ix) - 9^A"(x))]. 

Now we can consider a and b as free parameters that can be fixed from the requirement 
of the gauge invariance of the effective Lagrangian described by the one-nucleon loop 
diagram in Fig. 2c. 

Due to the constraints d'^Dj^{x) = d'^D^{x) = the corresponding terms in the 
Lagrangian (|3.14|) can be dropped 

/:Fig.2c(x) = te^[-bDiix)Dfsix)d''Af'ix) - 

~{b - a) Di{x) Dp{x) A"(x) - (3.16) 
-(a - b)d'^ Dl{x)D''{x)A(^{x) + b D"{x) d'^ D^{x) A^{x) + 
+ 3Dl{x)Dp{x){d''A^{x) - df^A^ix))]. 

The subsequent transformations are performed by applying the identity 

Dl{x)Di3{x)id''A^{x) - df^A^ix)) = (3.17) 
= df' Di{x)Dp{x)A''{x) - Dl{x)d''Df,{x)Af'{x), 

where we have dropped the total divergence and the terms proportional to 9" Dl^{x) and 
df^Dpix). 

Then it is convenient to rewrite the Lagrangian ( |3.16 ) as follows 



2 
V 

6 7r^ 

a - b)do, dUx) D"(a;) A^(x) + bDl{x) d" D^{x) A^{x) - (3.18) 



C-YiM^) = te^[-ia - b)Dl^ix)A^ix)D''{x) + b D^^ix) A(,{x) DUx) + 



-bDl{x)Dp{x)d" A^{x) - (6 - a)Di{x)Df^{x)d^ A'^ix 
+ 3Dl{x)Dp{x){d'' Af^ix) - d'^A'^ix))]. 
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Putting a = 26 we bring up the effective Lagrangian ( 3.18| ) to the irreducible form that 



contains two parts defining the renormahzation of the electric charge of the deuteron and 
the gauge invariant interaction coinciding with that given by Corben and Schwinger [14] 

/:Fig.2c(x) = ie-^^[-hDl^{x)A^{x)D^{x) + hD^'^{x)Ap{x)Dl{x) + 

+ (26 + 3) Di{x) Dp{x) F"^(x)]. (3.19) 

Now the effective Lagrangian CYi^.2c{x) is represented in the irreducible form, and we can 
proceed to the analysis of gauge invariance. 

Putting 6 = we remove the finite contributions to the renormahzation constant of 
the electic charge of the deuteron coming from the one-nucleon loop diagram in Fig. 2c 
and get a gauge invariant interaction. As a result the effective gauge invariant interaction 
reads 

Ccs{x) = ^ e ^ Dl{x) D,{x) F^-{x) . (3.20) 

Thus the one-nucleon loop diagram in Fig. 2c defines the effective Corben-Schwinger 
Lagrangian Cqs{x) describing the interaction of the deuteron with the electromagnetic 
field. 

We should underscore that the ambiguities connected with the shift of the virtual 
momentum of the diagram in Fig. 2c are the intrinsic peculiarity of fermion loop diagrams, 
caused by the computation of such diagrams within the cut-off regularization [15]. So, 
one cannot consider the parameters a and h connected to a shift of virtual momentum 
Q = aki + bk2 as parameters especially introduced in the model. The requirement of 
gauge invariance defines unambiguously the parameters a and b and, correspondingly, 
the effective Lagrangian Ccs{x). The choice a = 2 6 is unique, since it brings up the 
Lagrangian ( p.l8| ) to the irreducible form. The next choice 6 = is also unique, for the 
effective Lagrangian ( |3.19| ) is given by the irreducible form. 



One can show that the contributions of the anomalous magnetic dipole moments of 
the proton and neutron to the effective Corben-Schwinger Lagrangian are proportional to 
J2(Mn) and due to inequality Ad ^ Mn are small compared with that given by Eq. ( |3.20| ). 

4 The effective Aronson interaction in the relativistic 
field theory model of the deuteron 

In this Section we give the detailed derivation of the effective Aronson Lagrangian 
defined by the one-nucleon loop diagram in Fig. 2d and describing the effective electro- 
magnetic interactions of the deuteron. The effective Lagrangian described by the diagram 
in Fig. 2d is given by 

J d^xCFig.2dix) = j (fx j ^ ^ {2n)^^ D^p{x) Dl^{xi) ^(xa) x 

4 TT -^^^Y) 
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and 



1 



(4.2) 



xtr < a 



a/3 



k-Q M^-k-Q-h 



■7 



k — Q — ki — k2 



where the 4-vector Q = aki + bk2 is due to the shift of the virtual momentum of the 
nucleons in the diagram in Fig. 2d. The parameters a and b are free. They are not 
connected with those having been intriduced for the computation of the diagram in Fig. 2c. 

In the r.h.s. of Eq.( [4.1|) we have used the fields of the physical deuteron, for the 
renormalization introduces divergent terms that due to inequality Ad <^ Mn are 

small compared with the convergent ones we are keeping for the computation of the 
Aronson Lagrangian. The Gertsein-Jackiw ambiguity is given by 



6J''^^''\k,,k2;Q) = -tr(a"^Qa'^V)- 

o 



(4.3) 



Now we should proceed to the computation of J'°'^'^''^{ki, k2;Q). By analogy with ( |3.2| ) 
we get 



J'''''"'\k^,k2;Q) 



d'^k 



a(3 



M^ + k + Q 



M^-k^ 



2k-Q 

M^f-k^ 



a 



M^ + k + Q + ki 
M2 -F 



X 



2k-{Q + ki 



Ml^-k^ 



1 



^M^ + k + Q + ki + k2 



M^-k"^ 



2k-{Q + ki + k2) 



Ml^-k^ 



d^k 



TtH (M2 - P)3 



+a"^a^V(A; + Q + ki + k2)] + a'^f'ik + Q)a^''{k + Q + k^)^^ x 



tr [a"^(A; + Q)a'"'-f^ + a'"'af"'{k + Q + ki)-f^ + 



a/3/ 



x{k + Q + ki + k2)} 1 + 

_ rd^k 1 
~ J 7r2 i {Ml - k'^f 

+a''^af"'^\Q + h + k2 
+ 2 



2 A; ■ (3Q + 2ki + ^2 



al3 nvi 



+CT°''ff'"'(3 + 2 *i + kih^ + (t"»(t'"-7*(3 Q + 2ki + k-i)] - 



a/B nu X/ 



A;V"^(3Q + 2A;i + A;2)a^V} = J^l^'^^h, k2; Q) + J^2)'''^ 



{h,k2]Q) 



Integrating over k we obtain 



(4.4) 



J(ir'(A;i,fc2;Q) = -[l + 2J2(MN)]tr(a"^Qa'^V) + 
+ 1 tr K^a'^^(Q + fci)7^ + a^''a^''j\Q + h + h)] , 



(4.5) 
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^ tr K'^(3g + 2h + A;2)a^'^7^ + (4-6) 
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+a"^a'^"(3g + 2h + A;2)7 + fT"^a^V(3g + 2A;i + ^2)]. 
Now we should summarize the contributions and collect like terms 

J''^^'\hM\Q) = ^trK'5(g-i(2A;i + fc2))a^V] + 

o o 

+ i tr \a''^a^"'{kl - ^2)7^] + - tr [(T"^a'^^7^(A;i + 2 ^2)] 
6 6 



(4.7) 



It is seen that the g-dependence coincides with that obtained by the Gertsein-Jackiw 
method. Due to the arbitrariness of Q the vector {2ki + /c2)/6 can removed by the 
redefinition of Q. Thereby we get 



J''''^''\h,k2;Q) = -tr(a"'3ga'^^7^) + 





+ 1 tr a'^'^ik - h)!""] + 7 tr K^a^^7\A;i + 2 ^2)]- 
o 



After the calculation of the traces of the Dirac matrices we obtain J'"'^'^'^^{ki, k2', Q) leading 
to the following effective Lagrangian 



'^FiK.2d(a;) 



-le) 



9\ 1 



47r2 Ml 



o o 
+ ^ ih + a) Dlix) D^\x) Ax{x) + ^ (6 - a) DI{x) D''\x) d, A'^{x) - (4.9) 
16 



g Dlix) D''\x) Ax{x) + 8 Dl^{x) D''\x) (9^ Ax{x) - dx A^(x)) 



For the derivation of the effective Lagrangian (|4.9|) we have used the equation of motion 

dxD^^ix) + d^D,x{x) + d.Dx^^ix) = 0. (4.10) 

The analogous equation of motion is valid for the conjugated field. By collecting like 
terms in (|4.9|) we get 



le) 



9\ 1 



47r2 Ml 



{b + a - 1) DUx) D-'^x) Ax{x) + 



+ -{b - a - 3)DUx)D''\x)dxA^'{x) + 



(4.11) 



+-adxD^^''{x)D,^{x)A^{x) + - a Z}tM-(x) 5a I5^^(x) A^(x) 



The third and the fourth terms can be reduced by applying the equation of motion 

dxD^^ix) = -MlD^^x) (4.12) 
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and analoguous for the conjugated field. Then putting b + a 
obtain b 
form 



6 + a + 3, 



we 



2, that brings up the effective Lagrangian ( |4.11| ) to the following irreducible 



'^Fig.2d(a^) 

2 0t 1 . 



3 7r2 Ml 



Dl{x)A^{x)D''{x) - D^\x) A,{x) Dl{x) 
a) Dl.ix) D^^ix) Fx^ix). 



+ 



(4.13) 



Putting a = we remove the finite contributions to the renormalization constant of the 
electric charge of the deuteron and gain the gauge invariant interaction 



^Fis.2d(a^) 



z e 



2gl 1 



37r2 Ml 



Dl{x)D^\x)F; 



(4.14) 



Due to the irreducibility of the Lagrangian ( |4.13| ) the choice a = is unique. 

The effective Lagrangian ( [4. 141) coincides fully with that suggested by Aronson [10]. 
This means that the one-nucleon loop diagram in Fig. 2d defines unambiguously the effec- 
tive Lagrangian if one imposes the requirement of the gauge invariance. The same result 
can be gained if one requires the vanishing of finite contributions to the renormalization 
constant of the electric charge of the deuteron. 

By analogy with the computation of the Lagrangian ( [4. 141 ) one can obtain the con- 
tribution of the anomalous magnetic dipole moments of the proton and neutron to the 
effective Aronson interaction 



z e 




47r2 Ml 



DUx)D^\x)F,^{x). (4.15) 



As a result the complete expression of the Aronson effective Lagrangian is given by 



I e 



3 7r2 
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Ml 



DUx)D^\x)F^^{x). (4.16) 



Thus we have shown that in the relativistic field theory model of the deuteron and in one- 
nucleon loop approximation one can unambiguously compute effective electromagnetic 
interactions of the deuteron in terms of anomalous magnetic and electric quadrupole 
moments. This can be obtained under requirement of gauge invariance to every one- 
nucleon loop diagram separately. 

5 Radiative neutron— proton capture 

At low energies, the process of the radiative neutron-proton capture n + p — >■ D + 7 
comes off via electric and magnetic dipole transitions. In the usual notations ^^^^Lj, the 
deuteron is a 1 state, and the possible transitions are 



'Si ^ ^'Si (/i) 



0,1,2 



^Si(rf) , iSo^='Si(/i) 
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where /i = magnetic dipole and d = electric dipole. If the energies arc low enough, 
the nucleons are in the S-wave state and thereby only magnetic dipole transitions are 
possible. 

In the S-wave state the magnetic dipole moment operator acts only on spin variables. 
This implies that the transition ^Si ^Si is forbidden [16]. Thereby the only allowed 
transition is ^Sq ^Si. This means that the anomalous magnetic dipole moments of the 
proton and neutron should give the main contribution [16]. 

In the relativistic field theory model of the deuteron the interactions of the deuteron 
with other fields should be obtained through one-nucleon loop diagrams. Therefore, for 
the computation of the effective Lagrangian of the radiative neutron-proton capture we 
have first to compute the effective Lagrangian, describing low-energy neutron-proton 
scattering. Keeping to the one-pion exchange we obtain [1] 



= j^{[Pi^)n%x)][n<^{x)p{x)] + [pixh'n^ix)][n-{xh'p{x)] 



4M2 

+\p{x)^fj_^^n''{x)][n''{x)^''^^p{x)] + 3\p{x)^ i^n''{x)][n''{x)y'^p{x)] + (5.1) 
+^\p(x)a,,n'^{x)][n'^{x)a^''p(x)] }. 

Only terms [p{x)'~f^n'^{x)][n'^{x)'y^p{x)] and [p{x)'~ff^'~f^n'^{x)][n'^[x)^'^'y^p{x)] contribute to 
the S-wave of the neutron-proton scattering in the low-energy limit. Therefore, these 
terms should effect the ^Sq — >^ ^Si transition in the radiative neutron-proton capture at 
low energies. The corresponding one-nucleon loop diagrams are depicted in Figs. 3 and 
4. 

First let us consider the contribution of the diagram in Fig. 3a. The corresponding 
Lagrangian reads 

J d^X CFig,3^{x) = 

^ „-iki-xi -ik2-X2 ^iiki+k2}-x ■ ^ 9ttNN 9y_ n■^lv/, i . ^\ 

xe e e 327r^ \'^i-> W ) ■> 



where 



jr{k,M) = (5.3) 

d^k { . 1 ^ 1 ^ 1 1 

— tr<^75 . . . 

T^^i [ M^-k-Q M^-k-Q-ki - k - Q - ki - k2 ) 

and Q — a /ci + 6 /c2 is an arbitrary shift of virtual momentum. Fortunately, the integral 
J^^iki-i k2] Q) does not depend on the shift of the virtual momentum and can be computed 
unambiguously 

jriKk^.Q) = ^e'^'^'^^k.^k^p {e^^^^'^l), (5.4) 
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where we have used the relation tr(7^7^7'^7°7^) = —Aie'^'^"'^. The structure function (5.4) 
leads to the effective Lagrangian 

^Fig.3a(x) = - ^ DlixTF'^^ix) [n^{xh'p{x)] , (5.5) 

where *F^"'{x) = ^e'""'/^ F^p{x) . 

Now let us consider the contribution of the diagram in Fig. 3b. This contribution is 
caused by the anomalous mangnetic moment of the proton. The effective Lagranigan is 
defined 



J c/^x£Fig.3b(a;) 



/ ^-^^ Mx)j'p{x)] Dlix,) F^,{x,) X (5.6) 

2 



where 



jr^{h,h:Q) = (5.7) 

— tr <^ 7^ . ^ 7'^ . . ^ a"^ ^ . . ^ . 

Unfortunately, the integral J^'^^{ki^ k2; Q) depends on the shift of the virtual momentum. 
Therefore, the leading contribution in the momentum expansion is fully arbitrary and is 
given 

jr\kr,k2,Q) = -2te^''''''Q, = -2z£^"^^(afci + bk^)., (5.8) 

where a and b are arbitrary parameters. The contribution of the momentum ^2 produces 
in the Lagrangian the operator d ^* F ^'^{x) that vanishes due to Maxwell's equation of 
motion, i.e. d p*F ^^{x) = 0. Thus only the contribution of the momentum ki matters. 
The effective Lagrangian produced by J'^'^^{ki, k^, Q) given by (|5.8|) reads 

/:p,g,3b(x) = ^ «:p ^ DlixYF^^ix) [n^{x)^'p{x)] . (5.9) 

The effective Lagrangian corresponding to the diagrams in Figs. 3a is given by 

/:F.g.3a,b(x) = - (l - i '^p) DlixTF^'^ix) [n^(x)7^p(x)] . (5.10) 

Now let us discuss how we can fix the parameter a. As has been mentioned above the 
radiative capture n + p^D + 7at low energies is a magnetic dipole transition ^Sq 
^Si (/x). This implies that the contribution of the proton and neutron to amplitude of 
the radiative capture n + p — >^ D + 7 should be proportional to their magnetic dipole 
moments. In the case of the proton the total magnetic dipole moment equals (1 + Kp), 
i.e. /ip = 1 + Kp, whereas the total magnetic dipole moment of the neutron is just 
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its anomalous magnetic dipole moment, i.e. = /in. The effective Lagrangian ( ^.10| ) 
describes the contribution of the proton, therefore the quantity (1 — aHp/2) should be 
nothing more than the total magnetic dipole moment of the proton fip. This fixes the 
arbitrariness of the contribution of the diagram in Fig. 3b. Thus, we have to put a = —2. 
As a result the complete contribution of the diagrams in Figs. 3a, b reads 

/:Fig.3a,b(x) = - /xp ^ DlixyP'^^ix) [n^{xh'p{x)] . (5.11) 

Thus we have fixed the ambiguity introduced by the diagram in Fig. 3b by applying the 
selection rule ^Sq — >■ ^Si (/i). 

By extending the procedure of the computation of the diagrams in Figs. 3a,b on other 
diagrams in Fig. 3 we get the following complete effective Lagrangian 

U,,,,{x) = - (-"P - -"n) ^ DlixYF^'^ix) [n^{xh'p{x)] . (5.12) 

The contributions of the diagrams describing the interaction of the deuteron with the 
anomalous magnetic dipole moments of the proton and neutron via the tensor nucleon 
current, i.e. proportional to the constant g^, are divergent and due to the inequality 
Ad <^ Mn are small compared with the convergent ones. Thereby the contributions 
proportional to the coupling constant qt can be neglected and do not appear in the 
amplitude of the radiative neutron-proton capture n + p — >^ D + 7. 

The contribution of the [p{x)'jfj^'j^n'^{x)][n'^{x)'-f'^'-f^p{x)] interaction can be computed 
by analogy to that given above. As a result we obtain 

£pig,4(x) = -^e{^^p-^^^)j^^^Dl{xrF'^'^{x)Mxh,^'p{x)]■ (5.13) 

The total effective low-energy Lagrangian, describing the radiative neutron-proton cap- 
ture, is defined by the sum of the effective Lagrangians ( [5.12[ ) and ( |5.13| ) and reads 

C.s{x) = - (/^P - /^n) ^ DlixYF^'^ix) [n^{xh'p{x)] + 

2 " (5.14) 
(/ip - -"n) ^ 1^ DlixYF^-^ix) [n^(x)7.7'p(x)]. 

The amplitude of the radiative neutron-proton capture is given by 

A^(n + P^D + 7) ^ -(.p-.n)^^^."^-^.eKfc)e;(Q)x ^^^^^^ 

x[u-{pi){2Q, - Mn7.)7'm(p2)], 

where e^(/c) and e*^{Q) are the polarization 4-vectors of the photon and the deuteron, 
respectively, then u'^{pi) and u{p2) are the bispinorial wave-functions of the neutron and 
proton normalized by u'^{pi)u''{pi) = — 2Mn and u{p2)u{p2) = 2Mn. 

The cross section a{n + p — D + 7) of the radiative neutron-proton capture has been 
measured for thermal neutrons at laboratory velocities v/c = 7.34 -10"^ (the absolute 
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value is f = 2.2 ■ 10^ cm/sec) [17] 

or(n + p ^ D + 7)exp = (334.2 ± 0.5) mb . 
The cross section cr(n + p ^ D + 7) is given by 

(T(n + p-.D + 7) = -— / |-M(n + p ^ D + 7)|2 x 



(5.16) 



x(27r)V^)(Pi+P2-Q-A;) 



(5.17) 



(27r)32E^ (277)32^0 



where Ei, E2, and Ed are the energies of the neutron, proton, photon and the deuteron, 
respectively, v is the laboratory velocity of the neutron, and |7W(n + p — * D + 7)^ is the 
squared amplitude averaged over polarizations of the neutron and proton and summed 
over polarizations of the deuteron and photon 



|X(n + p^D + 7)|2 



4M2 



fi'TrNN 



Ml 
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2567r4 
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X - X 
4 



xtr 



(pi - Mn)(2Q, - Mn7.)7'(P2 + Mn)(-2Q- - Mn7")7' 



(5.18) 
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X [4(pi ■ p2 + 2M^) Q,Q^ + M^(pi,p- + p>2.) - (Pi ■ P2 - M^) g^^] , 

where a = /An = 1/137 is the fine structure constant. 

In the low-energy limit when the 3-momenta of the neutron and proton tend to zero 
we obtain 



|A<(n + p ^ D + 7)|2 = (/ip -^n)- 



1287r 



MI 



^XPpuj kak X 



X [4(pi ■ p2 + 2m2) g,g- + mUpiuP2 + P1P2.)] 



(^p - ^n)' 



25a Qp 
327r 



Ml 



(5.19) 



where we have used the relation ( p. 16 ). 

Thus the cross section of the radiative neutron-proton capture reads 



c,(n + p^D + 7) = -(/xp-/x„) 



■ 2 n2 

M2 



(5.20) 



Putting V = 7.34 -10 ^ we get the following theoretical value of a{n + p — > D + 7) 

cj(n + p ^ D + 7) = 156.4 mb . (5.21) 



It is seen that this theoretical value of the cross section disagrees with the experimental 
data. 
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However, it should be stressed that we do not have taken into account the resonance 
contibution to the ^Sq low-energy neutron-proton scattering [18-20]. The account of the 
resonance contribution gives the amplitude of the radiative neutron-proton capture in the 
form [21] 



A^(n + p ^ D + 7) 



x{u-{pi)[2QJl 



IQti Ml as 



- Mn7.]7'«(P2)}, 



(5.22) 



where as = — 23.748fm [4] is the ^Sq neutron-proton scattering length that is fully due 
to the resonance contribution [18-20]. 

As a result the cross section of the radiative neutron-proton capture is given by 



cr(n + p 
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IGttMI as 
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9nNN 
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IQuMl as 



(5.23) 



Mn^d = 308.8 mb. 



The theoretical value agrees now with the experimental data within an accuracy better 
than 8%. The cross section (|5.23|) predicted in the relativistic field theory model of the 
deuteron agrees also well with that given by the potential model [22]. 



6 The low— energy p + p^D + e+ + z/e reaction 

The process of the two-proton fusion p + p— >D + e+ + i/e plays an important role 
for the synthesis of deuterons in stars, where the deuterons are destroyed again by the 
reaction p + D ^ ^He + 7 [23]. In nuclear physics the process p + p — D + e"*" + z/e is 
a Gamow-Teller transition governed by the weak axial-vector nucleon current [24]. 

In the relativistic field theory model of the deuteron the low-energy process p + p — > 
D + e"*" + z/e is closely connected with low-energy proton-proton scattering, i.e. p + p 
— i> p + p. Low-energy proton-proton scattering differs very much from the low-energy 
neutron-proton scattering. This is mostly due to the strong contribution of the Coulomb 
repulsion [25]. The small relative velocities v in the proton-proton system are suppressed 
by the factor exp(— 27ra/f ), where a = e'^/An = 1/137 is the fine structure constant. The 
factor exp(— 27ra/f) is the Gamow penetration factor [25,26]. For the acquaintance with 
the earliest computation of the cross section of the p + p^D + e"^ + z/e scattering we 
refer readers to the paper by Bethe and Critchfield [27] and the book by Rosenfeld [28]. 

For the computation of the effective Lagrangian describing the low-energy proton- 
proton scattering we assume the one-pion-exchange approximation 

^l^i^) = - %# iPixh'pix)] [p{xh'p{x)] . (6.1) 



Here we have denoted gnNN^v) = gnNNC{v) and C{v) = y 27ra/t> exp(— yra/f) takes into 
account the Coulomb repulsion between protons at low energies [27,29]. 
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Applying the Fierz transformation we bring up the Lagrangian (|6.1|) to the form 



,2 



^cff(^) = {[p{x)p%x)mx)p{x)] + [p{x)jY{x)W{xh'pix)] + 

+ \p{x)-ff,-fY{x)W{x)-f''j''p{x)] + ^\p{x)a^,p%x)]\p'^{x)a^'''p{x)] }. (6.2) 

The process p + p— s^D + e^ + z/o should run via the intermediate W-boson exchange, 
i.e. p + p D + W+ — > D + e"*" + z/g. The Lagrangian describing the electroweak 
interactions of the W-boson with proton, neutron, positron and neutrino reads [2,19] 

^Z{x) = -0=[p{x)r{l - 9Al'Mx) + iye{x)Y{l - l')e{x)] W,{x) + h.c. . (6.3) 

Here is the electroweak coupling constant connected with the Fermi constant Gp = 
1.166 X 10-5 GeV-2 and the W-b oson mass Mw by the relation [2,19] 



F 



8M^ 72 



(6.4) 



where qa = 1.260 ± 0.012 is the axial-vector coupling constant [4] describing the renor- 
malization of the weak axial-vector hadron current by strong interactions. 

Only the interaction [p{x)'yn'j^p'^{x)][p'^{x)''y'^'y^p{x)] gives in our approach the main 
contribution to the amplitude of the transition p + p — > D + W. The corresponding 
one-nucleon loop diagrams are depicted in Fig. 5. 

The effective Lagrangian defined by the diagrams in Fig. 5 is given by 



J d'^x Cpig.^ix) 



/ ^-J0^ mxh.l'pix)] Dlix,) W^{x,) X (6.5) 



where 



J'""'ih,h;Q) = (6.6) 
^tr ^-^Y ^ ^ . ^ 7-7^ ^ ^ 



M^-k-Q M^-k-Q-ki - k - Q - ki - k2 



The computation of this integral can be reduced to the computation of the integral (|3T 
defining the Corben-Schwinger Lagrangian, i.e. 

J'^^'ik,, k2; Q) = J^^'ik,, k2; Q) + :r"^"^(A;i, k^, Q) , (6.7) 

where we have denoted 

J'^^'ik.M-.Q) = (6.8) 

J TxH M^-k-Q^ M^-k-Q-h^ M^-{k + Q + ki + fcs)^ J ' 
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The integral J^"^'^{ki,k2',Q) does not depend on the shift of virtual momentum, i.e. on 
Q, and in leading order in the long-wavelength expansion is given by 



tr {7^^ ^ 



ttH M^-{k + Q + ki + k2)^ 
1 .,1 



X 



^ ^ ^ 

M^~k-Q M^-k-Q-k 



-7 



d^k Mk 



1 + 



2k ■ (3g + 2ki + fca) 
M2 -P 



X 



7r22 (M2 - P)3 

tr{7°(MN + A; + Q)7^(Mn + k + Q + = 

d^k 2A/f^ - ^ ^ 

^ ■tr{7°Q7'^7'^ + 7V(Q + A:i)7"} 



7r2i (M^-fc2)3- 

^ (M2 -fc2)4 W7"(3g + 2A;i + A;2)7'^7'^ + 7"7^(3Q + 2A;i + k^h"} 
= -tr{7"Q7V + 7"7''(Q + ^1)7'} + 

+tr{7"(Q + Ih + Ikh^l" + 7 V(Q + pi + ^2)7^^} = 
= ^tr{7"(2fci + A:2)7''7' + 7"7''(-A:i + kh"} = 



(6.9) 



4{k^g^''-k'(g"^) + -{k^ + 2k2rg' 



Summing up the contributions (|6.9| ) and (p.l4|) we obtain 

- 10 2 2 

J'-'^ik,, k,- Q) = -{k'^g'^^ - k^gn " 3 (^1 + ^^Td'"' - ^k^g^'' + 

+\[{Q -ki + k^Yg'^' + (g - A;i + k^rg""^ +{Q-k^ + k^Yg^ ■ 



(6.10) 



In order to fix ambiguities produced by the shift of the virtual momentum let us note that 
k2, being the momentum of the W-boson, coincides with the momentum of the leptonic 
pair and is small. Therefore, we suggest to keep only the terms leading in ki expansion 

- 10 2 2 

J'^^^ik,, k^; Q) = -{k'^g"" - klgn - -(2 - a) k'^g^^ + -{a - 1) k^g'^'^. (6.11) 

Here we have dropped the term proportional to k^ giving in the Lagrangian the total 
divergence of the deuteron field that vanishes: d'^Djj^ = 0. 

The structure function jf"^''(fci, A;2; Q), defined by the Eg. ( |6.11| ), would be invariant 
under gauge transformations of the deuteron field if a = 3/2, i.e. 



J''^''{k,,k2;Q) = ?,{k'^g''^'-k'ig 



(6.12) 



The effective Lagrangian defined by the structure function (|6.12|) has the following form 



^FigA^) = -gA 



9w gl^M Sg 
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(6.13) 
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In turn the effective Lagrangian describing the low-energy process p + p— s^D + e^ + i/g 
reads 



V 



y/2 Ml 647r2 



DUx) [p'=(x)7Vl'(a;)], 



(6.14) 



where = z/c(a^)7^t(l — 7^)e(x) is the leptonic electrowealc current. 

In the low-energy limit when the 3-momenta of protons tend to zero the effective 
Lagrangian (|6.14 ) can be reduced to the expression 



£eff(a;) =igAM^ 



V2 Ml 327r2 
where we have used the relations 

\p''{x)Yl^p{x)\ -> -5(''°[p'=(x)7^p(x)] 
Dl,{x)^-iM^Dl{x) 

that is valid in the limit of low 3-momenta of protons and the deuteron. 
The amplitude defined by the effective Lagrangian ( |6.15D reads 



(6.15) 



(6.16) 



A^(p + p ^ D + e+z/e 



>CM9aM,-^^3^x ^^^^^ 



where u{ku) and f (/Cg) are the Dirac bispinors of the neutrino and positron. We have 
separated the factor C(f), describing the Coulomb repulsion, that is convenient for the 
sequent calculations. 

The cross section of the low-energy p + p— ^D + e^ + t'e reaction is defined 



cr(p + p^D + e+ + z/e) = - / |-M(p + p ^ D + e+ + i/e)|2 x 



X {2ny 5^ (pi + p2 - g - fce - K) jy^^ J't 



(6.18) 



(27r)32ED (27r)32Ee (27r)32E^ 



where u is a relative velocity of protons, and |A^(p + p — » D + e+ + z/e)P is the squared 
amplitude averaged over polarizations of protons and summed over polarizations of final 
particles, i.e. 



|A1(p + p ^ D + e+i/JP = C\v) giM, 



2 T^ji ^GlQ-p 

N 010^2 



Ml 



xtiiCke - me)7"(l - 7')A;.7''(1 - 7')} x ^ x 
X 1 X tr{(pi - Mn)7'(1 + l'h){P2 + Mn)7'(1 + I'h)}, 



9af3 



QaQp 



X 



(6.19) 



where me = 0.511 MeV is the mass of positron, si and S2 are 4-vectors of polarization of 
interacting protons such as sl = s\ = —1 and si ■ pi = S2 ■ P2 = [2]. 
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Since protons interact at low energies, we calculate the traces over Dirac matrices in 
the low-energy limit 
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tT{h-f"k-f^} = -24 ( EeE, - ■ K 



(6.20) 



X tr{(]3i - Mn)7'(1 + 7'si)(p2 + Mn)7'(1 + l^'h)} = M^{1 - si ■ S2). 
16 2 



Substituting Eq. (|6.2q) in Eq. (|6.19|) we get 



27GpQy) [fi'^NN 



M2 



X 



(6.21) 



X ( E^E^ - -fee -K Ul-si- S2] 



Now let us take into account the spinorial properties of the interacting protons. At low 
energies the two-proton fusion proceeds through an intermediate S-wave state with the 
total spin of the two protons equal to zero (S=0). The former means that the spins of the 
interating protons are antiparallel, that gives Si-S2 = —1- However, Si-S2 = —1 can be also 
realized for the two-proton state with a total spin S=l and zero projection onto the z-axis. 
Therefore, the factor (1 — si-S2)/2 sould be expanded as (1— si-S2)/2 = (l/2)s=o+(l/2)s=i. 
Keeping only the S-wave contribution we arrive at the expression 



|M(p + p^D + e+z/e)|2, 
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(6.22) 



X I E^Ey — — fcg ■ fc^ 



Now we should carry out the integration over the phase volume of the final D e^i^e^state 
d!^Q d?K d^K 



(27r)32ED (27r)32Ee {2Tif2E^ 
X (27r)^ 5^^\pi +V2-Q-K-K) (e,e, - ■ K ^ 

1 /"^D j 

^ El - mlE,{e^ - E,f dE 



(6.23) 



327r3MN 7™, V ^ " " ''^ ' 9607r3MN 
where ^ = meje^ and the function /(^) is defined by the integral 

\J — S,"^ X (1 — x^dx = 



30 



9 



15 



0.776 



(6.24) 
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and normalized to unity at ^ = 0, i.e. /(O) = 1. 

The cross section of the S-wave low-energy two-proton fusion p + p^D + e^ + z/g 
is given by 



cT(p + p^D + e+ + i/e)s - 



Ml 



1.88 X 10-^9 ^!Mcm2. 



(6.25) 



The cross section is calculated in units of ^ = c = 1. The appearance of the factor C^(f ) 
agrees well with the result obtained by Bethe and Critchfield [27]. 

Now we should average the cross section over relative velocities of the proton-proton 
system with a Maxwell-Boltzmann distribution [27] 

< i; cr(p + p ^ D + e+ + z/e)s >= 

3/2 



1.88 X 10-^^ - { f Sv C\v) exp(-^v 



in-49, 2 ( MmY' / Mr, 2 27ra\ (6.26) 

1.88 x 10 4.«U^ dvve.pi-^v --- = 



where k = 8.62 x 10~^^ MeV ■ is the Boltzmann constant, and T is the temperature. 
Following [27] and computing the integral by a saddle point approximation we get 



(6.27) 



1 2 flV 

< w cr(p + p D + e+ + i^e)s >= 1-46 X 10^''° x - x - x - x 

a 71 \3 J 

X X r^e^^cm^s"^ = 1.45 x lO^^V^ e"^ cm^ s"\ 
where r is connected with the temperature [27] 

-3 ^ • 



The temperature dependence of Eq. ( |6.27| ) coincides with that obtained by Bethe and 
Critchfield [27]. 

Setting T = Tc = Tg = 15.5, measured in unites of 10^ K, where Tc is the temperature 
of the solar core in the Standard Solar model [24], that gives r = 13.56 we get the following 
estimate 

< (t(p + p ^ D + e+ + z/c)s > = 3.45 x 10"^^ cm=^ s"^ . (6.29) 

This magnitude is larger by a factor of 2.9 compared to the potential approach [24] (see 
also [30]). In order to reconcile our result with the solar luminisoty we should assume that 
the temperature in the solar core equals Tc = 13.8 in 10^ K [31]. As has been remarked in 
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Ref. [31] the enhancement of the magnitude of the cross section of the two-proton fusion 
leads to a decrease of the temperature in the solar core. This gives a strong suppression 
of the solar neutrino fluxes. 

The increase of the amplitude of the p + p — D + W transition found in our ap- 
proach is connected with the computation of the amplitude in terms of one-nucleon loop 
diagrams. Indeed, the structure function Eq. (|6.12 ) defining the effective Lagrangian of 
the p + p — s> D + W transition is due to the contribution of the anomalous part of the 
one-nucleon loop diagram of the AAV-kind, i.e. with two axial-vector and one vector 
vertices [15,32]. One cannot expect that such an anomalous contribution should coincide 
with the contribution described within potential approach in terms of the overlap integral 
of the wave-functions of the deuteron and two-protons. The ambiguity of the anoma- 
lous contribution has been fixed by the requirement of gauge invariance under the gauge 
transformations of the deuteron field. It is very similar to the fixing of the ambiguity of 
the well-known axial Adler-Bell-Jackiw anomaly [2,33]. 



7 Conclusion 

In the present paper we have developed the relativistic field theory model of the 
deuteron that has been suggested in Ref. [1]. We have given the elaborate computation 
of the Corben-Schwinger and Aronson effective Lagrangians describing the interactions of 
the deuteron with electromagnetic fields and defining the anomalous magnetic and electric 
quadrupole moments of the deuteron. We have adjusted the model to the calculation of 
the low-energy processes such as the radiative neutron-proton capture n -|- p — D + 7 
for thermal neutrons and the fusion of two protons p + p— >D-|-e+-|-z/e. This reaction 
is very important for stellar nucleosynthesis, where the deuterons are destroyed again by 
the reaction p + D — ^He -|- 7 [23]. Unfortunately, our model is far from being induced 
by the dynamics of QCD and seems like an old-fashion approach. It is based on the 
Lagrangian Eq. ( |2.1| ) written with the only guiding principle that the deuteron is a bound 
state of the proton and neutron with spin one [3,10]. 

We have shown that the model, supplemented by the dynamics of low-energy neutron- 
proton and proton-proton scattering, is able to describe reasonably the cross sections of 
the radiative neutron-proton capture in agreement with both experimental data and 
predictions of the potential model. For the cross section of p + p ^ D + e+ + z/e we 
have obtained an enhancement by a factor 2.9 compared to the potential approach. This 
enhancement of the cross section could be of relevance for a suppression of the solar 
neutrino fluxes [31]. 

We have found that the coupling constant does not contribute to the amplitudes 
of the processes under consideration. Thereby one can assume that the interaction of the 
deuteron with the antisymmetric tensor nucleon current is less important for the physics 
of low-energy interactions of the deuteron than the interaction with a vector nucleon 
current. If it is true, this gives rise the question concerning the need of keeping a nonzero 
value for the coupling constant g^- What would happen if we would put = 0? 

Putting (yfx = we face the only problem of the fit of the binding energy of the 
deuteron by the cut-off Ad = 45.688 MeV. Using Eq. (|2.18| ) at = we obtain the 



best fit of the binding energy at Ad = 64.843 MeV. Recall that we have identified 1/Ad 
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with the effective radius of the deuteron [3]. By using the new value of the cut-off we get: 
td = 1/Ad = 3.043 fm. This value agrees well the average value of the deuteron radius, 
i.e. < r > = 3.140 fm [34]. 

Putting then qt = we simplify our model reducing the number of free parame- 
ters and gain the cut-off Ad = 64.843 MeV that should define the average value of the 
deuteron radius td = 1/Ad = 3.043 fm. 

Thus the number of free parameters introduced for the definition of parameters of the 
physical deuteron are just three: gy, Ad and A and we could fix them unambiguously. In 
this connection we should underscore that additional free paprameters appearing due to 
the ambiguities connected with the shifts of virtual momenta in the diagrams in Figs. 2-5 
are intrinsic peculiarities of one-fermion loop diagrams, caused by the computation of 
such diagrams within cut-off regularization [15]. Therefore, one cannot consider these 
parameters as free parameters especially introduced in the model. Of course, the appear- 
ance of these ambiguities can lower, from the first glance, the predictive abilities of the 
model. However, as we have shown, most of these parameters can be unambiguously fixed 
by the requirement of gauge invariance. 

Finally, we want to estimate the theoretical uncertainty of the relativistic field theory 
model of the deuteron. Unfortunately, in our model there is not a small parameter, like 
a = e^/47r = 1/137, the fine structure constant in QED. There is not a large parameter 
like the number of quark colours N in the multi-colour extension of QCD with SU{3)c —>■ 
SU{N)c that allows to apply large A^-expansion in powers of 1 /N to quark-gluon diagrams 
describing the strong low-energy interactions of hadrons. 

Our approach is effective and based on the one-nucleon loop diagram approxima- 
tion for the description of both the self-interactions of the deuteron and coupling of the 
deuteron to other particles. By describing the self-interactions of the deuteron at the 
one-nucleon loop approximation we have fitted all parameters characterising the physical 
deuteron. In this case the only way to estimate the theoretical uncertainty of the ap- 
proach is to compare the one-nucleon loop data with the two-nucleon loop corrections. 
In Ref. [1] we have calculated the two-nucleon loop correction to the binding energy of 
the deuteron. Setting = 0, this correction reads [1] 

The numerical value of the two-nucleon loop correction makes up about 16% of the binding 
energy = 2.273 MeV calculated in the one-nucleon loop approximation (see Eq. ( ^.14D ). 
Thus the magnitude A = ±16% might be accepted as a theoretical uncertainty of the 
relativistic field theory model of the deuteron. This agrees well with our prediction for 
the cross section of the radiative neutron-proton capture. Indeed, the theoretical value 
(T(n + p — > D + 7) = 308.8 mb reproduces the experimental data a(ri + p — > D + 7)cxp = 
(334.2 ± 0.5) mb with an accuracy of about 8%. By taking into account the theoretical 
uncertainty the cross section of the radiative neutron-proton capture should read: a{n + 
p — s> D + 7) = (308.8 ± 49.4) mb. With regard to the two-proton fusion cross section, we 
predict an enhancement of 2.9 ± 0.5 compared to the potential approach. 

We acknowledge fruitful discussions with Prof. G. E. Rutkovsky and Dr. H. Leeb, and 
Dr. M. Meinhart for discussions concerning effective radius of the deuteron. This work 
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Figure Captions 



• Fig. 1 One-nucleon loop diagrams contributing to the binding energy of the physical 
deuteron, where rf — Cfi^ is the field of anti-neutron. 

• Fig. 2 One-nucleon loop diagrams describing the effective Corben-Schwinger and 
Aronson interactions that are responsible for the anomalous magnetic and electric 
quadrupole moments of the physical deuteron, where rf — CnF is the field of anti- 
neutron. 

• Fig. 3 The contribution of the \p{x)^^n'^{x)\[nF{x)^^p{x)\ to the amplitude of the 
radiative neutron-proton capture. 

• Fig. 4 The contribution of the [p(a;)7^7^n'^(a;)] [n'=(a;)7^7^]9(a;)] to the amplitude of 
the radiative neutron-proton capture. 

• Fig. 5 The contribution of the [p(a;)7^7^p'^(a;)][p'=(a;)7'^7^p(a;)] to the amplitude of 
the p + p— >D + e+ + i/e scattering. 
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Erratum 

to the paper "On the relativistic field theory model of 
the deuteron II" 

by Ivanov et al. published in Nucl. Phys. A617 (1997) 414. 



Abstract 

We correct the value of the cross section for pp-fusion (p + p ^ D + e"*" + i^e) 
calculated in Ref. [1]. We find a contribution to the astrophysical factor 6Spp{0) = 
2.01 X 10~^^ (1 lb 0.30) MeVb which is obtained only due to weak and strong low- 
energy interactions of the protons and the deuteron when neglecting the Coulomb 
repulsion between protons in the intermediate state. Minor misprints in Ref. [1] are 
corrected. 



1. The enhancement factor 2.9 obtained in Ref. [1] with respect to the potential model 
is mainly related to the neglect of the electromagnetic mass difference of the proton and 
the neutron by using a value of W = 2.223 MeV instead of 0.932 MeV in the phase-space 
factor {W + E)^ f {me/ {W + E)) . Furthermore, we used a few erroneous combinatorial and 
dynamical factors connected with identical and dynamical properties of the interacting 
protons. All of these factors are included now correctly in the Eqs. 1 and 2 given below. 
For the astrophysical factor we give the following new expression 



6Spp{E) 



a 



9glGlQr,M^ 
25607r4 



4M2 



8x/27rM^ 

i/TrNN 



"PP 
Mm 



{W + Effi 



W + E 



,(7.30) 



where E is the kinetic energy of the relative movement of the protons, and the function 
fimJ{W + E)) is given by Eq. (6.24) of Ref. [1]. For 5^pp(0) we obtain 



5S'pp(0) = a 



9glGlQj,M, 



25607r4 



4M2 



1 - 



8V2ttM^ 

i/TrNN 



Mm 




2.01 X 10"^^ (1 ± 0.30)MeVb, 



(7.31) 



where 30% is the assumed theoretical uncertainty of our model. The real theoretical 
uncertainty of the approach can turn out to be much less. The value 5S'pp(0) = 2.01 x 
10-25 (1 ± 0.30)MeVb makes up (51.7 ± 15.5)% of the value Sl^{^) = 3.89 x lO'^^ (1 ± 
0.011)MeVb obtained by Kamionkowski and Bahcall in the potential approach [2]. 

We stress that the astrophysical factor 6Spp in Eqs. ( [7.30| ) and (|7.31|) describes only 
the contribution of the weak and strong low-energy interactions and neglects the Coulomb 
repulsion between protons in the one-nucleon loop. The Coulomb repulsion in the initial 
state is included in terms of the Gamow penetration factor C{v) = ^/2^^] exp{—7iri) [3], 
where t] = a/v, multiplied by the amplitude defined by weak and strong low-energy 
interactions of the protons and the deuteron. 

Our result concerning the astrophysical factor 5S'pp(0) = 2.01 x 10^^^ (l±0.30)MeVb is 
rather promising. Indeed, as has been shown by Kamionkowski and Bahcall an important 



36 



contribution to the astrophysical factor for pp-fusion comes from the region where the 
Coulomb repulsion between the protons is dominating [2]. Since we have managed to 
describe about 50 % of the astrophysical factor obtained in the potential approach keeping 
only the contribution of weak and strong low-energy interactions of the protons and the 
deuteron, we expect to get the contribution of the same order when including the Coulomb 
repulsion in the one-nucleon loop. These calculations are in progress now. We note that 
we have tested our approach by calculating the disintegration of the deuteron by reactor 
antineutrinos Pe + D— >n + n + e+. This process is governed by the same physics as the 
quantity SSpp{0). We find excellent agreement with the potential model results and the 
experimental data [4]. 



2. Eq. (5.22) of Ref. [1] should read 

A<(n + p-^D + 7) = -in, - Hn) ^^j^ ^ e"""' K 4{k) e^Q) x 

x{«^(p0(2g.- MN7.)7VP2)}fl - 



Eq. (5.23) of Ref. [1] should read 



a(n + p^D + 7) = - (/ip - /^n) 



X 1 



N, 



5'7rNN 

M2 



Mn£d = (276±83)mb, 



X 



where = as in notations of Ref. [1]. Our result agrees reasonably well with the 
experimental data (T(n + p ^ D + 7)cxp = (334.2 ± 0.5) mb [1]. 

The change of the cross section for radiative captiu^c is caused by the equal contribu- 
tions of the pole on the unphysical sheet to the amplitudes defined by the 7^ x 7^ and 
7^7^ X 7/i7^ neutron-proton interactions at low energies. 

In Ref. [1] we have predicted an accuracy of the approach of about 16%. Since this 
accuracy has been obtained from the comparison of the one-nucleon loop and two-nucleon 
loop contributions to the binding energy of the deuteron, it should be understood as the 
accuracy of amplitudes. Thus,the theoretical accuracy of cross sections should be about 
30%. However, in practice this accuracy can be much better. 



3. Minor misprints: 

At the end of Chapter 2 of Ref. [1] the phrase "The contributions of these diagrams 
are divergent and can be removed by the renormalization of the wave-function of the 
deuteron" should read "The contributions of these diagrams are divergent and can be 
removed by the renormalization of the electric charge of the deuteron" . 

The phrase "... finite contributions to the renormalization constant of the wave- 
function of the deuteron ..." appearing above Eqs.(3.20), (4.14) and (4.15) should read 
"... finite contributions to the renormalization constant of the electric charge of the 
deuteron ..." 

In the r.h.s of Eq. (3.4) the common sign instead of (— ) should be (-I-). 
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In the r.h.s of Eq.(5.18) g'^ should replaced by M^g'^. 

Below Eq. (6.1) g^NNiv) = g^NNC{v) should read S'^nnI^) = 91nnC{v). 

In the r.h.s. of Eq. (6.23) tt^ should be replaced by tt^. 

Below Eq. (6.29) the phrase " This magnitude is smaller ... " should read " This 
magnitude is larger ... ". 
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